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Abstract 

Determination of Higgs self-interactions through the double Higgs production from gluon fusion 
is a major goal of current and future collider experiments. We point out this channel could help 
disentangle and resolve the nature of ultraviolet contributions to Higgs couplings to two gluons. 
Analytic properties of the double Higgs amplitudes near kinematic threshold are used to study 
features resulting from scalar and fermionic loop particles mediating the interaction. Focusing on 
the hh invariant mass spectrum, we consider the effect from anomalous top and bottom Yukawa 
couplings, as well as from scalar and fermionic loop particles. In particular, the spectrum at high 
hh invariant mass is sensitive to the spin of the particles in the loop. 
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I. INTRODUCTION 


Now that the Higgs boson has been discovered at rrih = 125 GeV, the next important 
task is a detailed exploration of the Higgs properties. The measured Higgs boson production 
rates and the extracted values of the Higgs couplings are close to the Standard Model (SM) 
predictions, but at the 0(10 — 20)% level, there is room for new physics effects in the 
Higgs sector. The structure of the Higgs potential is completely determined in the SM and 
measuring the Higgs self-interactions is an important step in determining if the observed 
boson is identical to the Higgs boson predicted by the SM. The Higgs self-interactions are 
most directly probed by double Higgs (2h) production at the Large Hadron Collider (LHC), 
gg —?■ hh, which has a very small rate, a 34 fb at = 13 TeV [M], making this 
measurement only feasible at high luminosity [81-1T8]. In the SM the cross section receives 
contributions from both box and triangle diagrams, and the large cancellation between the 
diagrams at threshold makes the gg —)■ hh process particularly sensitive to new physics 
contributions [TUH3U] . 

Beyond-the-SM (BSM) physics can contribute to 2h production in a variety of ways, 
including anomalous tth [bbh) and tthh [bbhh) couplings [311138] . resonant enhancements 
[391148] , exotic decays [HI |19] and new colored scalar [501155] or fermonic ^6H62] particles 
contributing to the loop amplitudes. Some of these effects, for example the modihed tth 
couplings or new colored particles in the loop, also affect the single Higgs (Ih) production 
in the gg ^ h channel. However, it is difficult to disentangle new physics effects in Ih 
production because of the limited number of kinematic observables in the hnal state. Using 
the higher-order process oi gg ^ h + j may help with measuring the top Yukawa coupling 
[63l l64] . but is unlikely to resolve the nature of the colored particle mediating the loop. In 
this work we are interested in the question of whether 2h production is sensitive to the un¬ 
derlying ultraviolet source of new physics and can potentially differentiate between different 
sources of new physics. We will see that many of the aforementioned new physics effects 
can signihcantly change the rate as well as the kinematic distributions in 2h production. 
(See Refs. |2Hl |52] for previous studies of new physics effects in the 2h kinematic distribu¬ 
tions.) In some cases, the changes are severely restricted by the (close to SM predicted) 
measurements of Ih production. 

This work is organized as follows. In Section we review the basics of Ih and 2h 
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production to set our notation. One of our major new results is in Section where we 
discuss the analytic structure of the 2 h amplitude near threshold in the case where the new 
physics arises from heavy fermions or from heavy colored scalars in the loops. Section IV 
contains numerical results for 2h production at y/S = 13 and 100 TeV. The amplitude for 
2 h production from intermediate colored scalar loops is reviewed in an appendix. 


II. BASICS OF Ih AND 2h CROSS-SECTIONS 

In this section, we review the lowest order results for Ih and 2h production from gluon 
fusion in order to £x our notation. We begin by presenting an effective Lagrangian, and 
then consider the specihc contributions from heavy fermions and heavy colored scalars. 

A. Non-SM Interactions 


We consider the following effective Lagrangian, where we are only interested in new 
physics affecting Higgs rates, assume SM kinetic energy terms {Lke)i and assume no light 
particles other than those of the SM. Including only third generation fermions. 




Leff = LkE ~ ~ 

_ f 1 + f 1 + 54 I - ^^G^’^"G7h2, (1) 

\ J 2v \ ) 8 v-^ 12 'kv 2A'kv^ 

where in the SM, 6 t = 6 b = Cg = Cgg = c^l ~ ^fh = 63 = 64 = 0. Global fits to Higgs 
production rates at the LHC limit the deviations of 6 t and Cg from 0 in a correlated fashion. 


as described below in Eq. Deviations of the 6-Yukawa coupling from the SM prediction, 
6 b, are less constrained 


B. Colored scalars 


The contributions from colored scalars depend on the parameters of the scalar potential. 
We use the following Lagrangian for an SU{2)l singlet, SU{3)c complex scalar, s. 


A. 


Cs,c = {Dgs)*{D^s) — mls*s — ^{3*3^ — ks*s\H^H 


—)■ {Dg3)*{D'^3) — ml3*3 — —{ 3 * 3 )“^ — K,3*3 


{h + v) 


72 


( 2 ) 

( 3 ) 
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where H is the SM SU{2)l doublet with {H) = (0,u/v^)^- In this normalization, the Fermi 
constant Gp = l/(v^n^) and v ~ 246 GeV. If the scalar, s, is real, 

= \(D,s)(D'‘s) - - F/ - \H'H\ ^ (4) 

The physical mass for either a real or complex scalar is, 

ms = mo + — , (5) 

where mo = 0 is the limit where the scalar gets all of its mass from electroweak symmetry 
breaking. The cubic and quartic scalar couplings are, 

r * 7 Q2h‘2s 7 2 

L ~ -gh2sS sh - —s sh 

gh2s = KV , g2h2s = K , ( 6 ) 

and similarly for a real scalar. 


C. Ih Production from Scalars and Fermions 


The leading-order (LO) gg —)■ Ih production rates due to virtual scalars and fermions are 
well-known. It is convenient to introduce the loop functions 


Fi/2ir) = - 2 r[l + (1 - r)/(r)] 
Fo(r) = r[l - r/(r)] , 


(7) 

( 8 ) 


where r* = 4:m‘^/m\ and 


/M = 


n 2 


sin-'(l/y?) 


if r > 1 


,1-1- \J\ —T 
In I - , I — n\ 


n 2 


if r < 1 . 


(9) 


,1 - yi-T) 

Then including colored scalars and non-SM fermion interactions as dehned in the previous 
subsections, 

-2 I 2(1 + ^;,) 


^ A) = 


a: 


10247r 


Zw.)- 


h 


Ti/7(r7.) + c,(-^) 




gh2siV „ , ^ 

mi. 


( 10 ) 
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FIG. 1: Top row: Triangle diagrams for gg —)• hh. Bubble diagrams with quartic scalar-gluon 
vertices are included with the triangle diagrams in our results. Bottom row: Box diagrams for 
gg ^ hh. Triangle diagrams with quartic scalar-gluon vertices are included with the box diagrams. 


where T(-) is the Dynkin index for the corresponding representation under SU{3)c defined as 
TriT^T^) = 6 abT{-), and 6 r = 1/2 for real scalars and 1 for complex scalars. The Dynkin 
index is 1/2 for fundamental representations and 3 for adjoint representations, respectively. 
For SM fermions, T(/) = 1/2. 

Neglecting the 6-quark contribution and noting that F 1/2 and Fq are well approximated 
by their large mass limits, Fi/ 2 (rt —>■ 00 ) = —| and Fq{ts —)■ 00 ) —)■ — 

^ (Tjgg h) 

^ ~ (^{gg -t h) \sM 

^ 1 + 2(^<^t + Cg + . (11) 


Eq. 11 is the well-known result that Ih production has little discriminating power between 
6 t and Cg [UTHTO] . The coefficients of Cg and Cgg from heavy colored scalars can be found in 
Refs. [7TH75] . In the SM, gg —)■ Ih production receives significant QCD corrections beyond 
LO QCD. The NNLO contributions from arbitrary fermions m and scalars [751 [76] are 
significant. However, since we are typically concerned with ratios relative to the SM, we 
work at leading order. 


D. 2h Production for Fermions and Scalars 

The LO gg —?• 2h production rates from fermion and scalar loops can be found in Refs. [U 
[2| and Refs. [52l[53|, respectively. The LO partonic cross-section for g{pi)g{p 2 ) —t h{ki)h{k 2 ) 
is given by 

2ft) = (12) 
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where s = (pi -hp 2 )^, i = (Pi — ^i)^- In the above 1/2^ comes from averaging over the initial 
gluon helicities, 1/8^ from color averaging and 1/2! from the identical hnal state particles. 
The amplitude-squared can be written as 


\Mf = 




+ <53)|5^(1 + 


s — ml 


fi 


+U + 


3ml{l + 63) 


s — ml 


I mi. 

^ Hi ** 




d/d 


_ 2^2 ^ ^Fa(s, f, m\, ml) + ^(1 + 5/J^Fn(s, t, ml, m\) 


h 




h 


E r rpf \9h2si'^ 17 r 2 2 \ ‘^^99 

5RT[Si) Fu (s, t, m;,, m,. ^ 


mz 


5^(1 + ml,ml) + 


t, ml, mil 


m 


Si 


.(13) 


where — 1 comes from summing over the gluon color index. The form factors 
F/\{s,t,ml,ml), Fn{s,t,ml,ml), and Gn{s,t,ml,ml) resulting from the SM top quark 
are given in the appendix of Ref. [I]^. The Feynman diagrams for the case of scalar particles 
are shown in Fig. [T| and the corresponding form factors are given in Appendix]^. In the 
large mass limits, 


F^{s,i, ml,ml) 

-P- 

2 1 ^ . 

3 ^ oo) 

FA{s,i, ml, ml.) 

-)■ 

1 1 ^ / 
g -t --Fo{ts -)■ oo) 

Fa{s,i, ml,m\) 

-)■ 

2 1^. 

-- -t oo) 

Fu{s,i, ml,m\) 

-)■ 

1 1 ^ , , 
-g -t 

Gu{s,i, ml,m^) 

-)■ 

o(^). 

V m^ / 


(14) 


^ In the SM, including only the top quark contribution Eq. 13 differs from Eq. 13 in Ref. [T] as well as 
Eq. 4 in Ref. [5^, but agrees with Eq. 6 in Ref. [5], Eq. 5 in Ref. m and Eq. 4 in [55], after plugging 
in Gp = l/(-\/2'c^) and taking account differences in the normalization of the form factors employed. 

^ We disagree with the overall normalization of the corresponding expressions in Refs. jSS] . In addition, 
the first 2 arguments of the last D function in Eq. 16 of Ref. [S3] should be swapped. 
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In the large mass limit, only the the spin-0 contribution survives, 
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( 16 ) 


The contributions from the anomalous couplings in Eq. 15 are consistent with those in Refs. 


[2HI EH ESI |37|. Furthermore, the contributions to Eq. ^ which come from the triangle 
diagrams are related to the fermionic and scalar contributions to the 1-loop QCD {3 function 
via the Higgs low-energy theorems ITU, which can be used to systematically compute higher 
order QCD corrections to the triangle loops [73l [78] . 


III. ANALYTIC STRUCTURE 

A closed-form analytic expression for the 2h production amplitude at threshold, s = 4m^, 
may be obtained from the imaginary part of the amplitude, combined with a knowledge of 
the amplitude’s limiting behavior as the particle masses in the loops go to either zero or 
inhnity |7S|. Alternatively, the threshold result can be obtained by a direct expansion of 
the full amplitude. For a colored fermion of mass mf running in the loops, the separate 
components of the amplitude arising from the triangle and box diagrams are, at threshold, 

\th = F^l\s = 4ml,i = -ml,ml,m‘j) 

\th = F^^\s = Aml,i = -ml,ml,m)) 

(16) 


where Tf = Am^/ml and T{f) is again the Dynkin index of the SU{3) representation of the 
fermion. The total amplitude is proportional to the sum of the two expressions above. 


= I - 1 + D (1 - Y) 


+ 


(? “ 0 ^ (^) ^ ~ - larcsin + 

(y - l) (U + 1) arcsinh^ “ (U “ 3)^7^Tlarcsinh 
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In the heavy mass limit, 


F'«U « -^r7>-*r7'^ + 0(ry) (18) 

and agrees with the result of ug. 

For a scalar of mass rUg with k = kq = 2ml/v‘^, the triangle and box amplitudes at thresh¬ 
old are found by analytic continuation of the imaginary contributions given in Appendix B, 



I K,=K,o 

\th 


T^{s) \K=Kq 

\th 


= Ami, i = -ml, ml, m)) 

1 

16 ^ 

^(■5) / g _ 


— T{s)ts ( 4 — Ts arcsin^ 
16 


Ts 

F^^>{s = Ami, i = -ml, ml, m)) 


(2 4- Ts) arcsin^ 
(2 -|- Ts) arcsinh^ 


+ 


Vf, 

+ 2 \/ts — 1 arcsin 
2 \/t 7^FT arcsinh 


+ 


(19) 


where = Ami/ml ^(■s) f^e Dynkin index of the scalar’s SU{3) representation. 

(s) 

We have included bubble diagrams with quartic scalar-gluon couplings in F)^ and triangle 

( s ) 

diagrams with such couplings in h Just as for fermions, we may expand the total threshold 
amplitude as 


rp { s ) |k=k;o |k=ko i p(®) 

^ \th \th \th 1 


( 20 ) 


where the functions F^'* \th^ \th are the threshold values of the form factors F^, F^^'^ in 


Eq. 15 In the heavy mass limit, the threshold result is 


F<”> irr” « -a' - + o(Tf ) 


( 21 ) 


The cancellations between the triangle and box functions for fermions and scalars are 
shown in Fig. In each panel, the cancellation clearly gets more exact for heavy loop 
particles. However, due to the small coefficients in the expansions above, the amplitude at 
threshold is still signihcantly suppressed for hnite masses. For the SM top, the indicated 
point in the left panel of Fig. shows that the triangle and box functions cancel to (9(10%). 
This cancellation will be spoiled by a non-SM Higgs self-coupling, additional interactions 
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FIG. 2: Cancellation between contributions to the 2h amplitude for fermions (left) or scalars (right) 
that get their mass entirely through couplings to the Higgs doublet. In the left panel, the value of 
Tf corresponding to the top quark is indicated. 



1+^3 


FIG. 3: 2h amplitude at threshold in the SM, as a function of the Higgs self-coupling. Only top 
quarks are included, with mt =173 Ge V. 

between the fermions and the Higgs boson, or if the scalar mass receives a contribntion that 
is not from the Higgs (he., k ^ kq). 

Having established that the cancellation between Fa \th and Fu \th is largely present 
for hnite loop particle masses, we now investigate the behavior of the cancellation in the 
presence of additional conplings. It is natnral to begin by considering the effect of a non- 
SM Higgs self-conpling on the 2h amplitnde from top loops at threshold. Snch a rescaling 
wonld affect \th only, since the box diagrams do not involve the Higgs self-conpling. 
Fig. shows how a modified Higgs tri-linear conpling wonld significantly alter the threshold 
amplitnde, leading to the well-known resnlt that 2 h prodnction is a sensitive probe of the 
Higgs self-conpling. Indeed, for arbitrary loop particle mass, there is a perfect cancellation 
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FIG. 4; 2h amplitude at threshold for scalars that get different contributions to their masses from 
electroweak symmetry breaking. The left panel shows curves of constant k/kq, while the right panel 
shows curves of constant physical scalar mass mg. 


between the one-loop triangle and box diagrams for 

1 + <^3 = 


p(/) I 

\th 


p(/) I 

■k A \th 


( 22 ) 


where the ratio of the box and triangle diagrams at threshold approaches -1 as the loop 
particle gets infinitely heavy. For the SM top mass, the cancellation is perfect when ^3 
0.09. The next section considers fnrther new conplings between the SM qnarks and the 
Higgs bosons. 

We now tnrn to the scalar case, inclnding arbitrary soft masses of scalars conpling to 
the Higgs. For a scalar which does not receive its mass entirely from electroweak symmetry 
breaking, mo 7 ^ 0 , the amplitnde at threshold is 


1 _ ^ p(s) |K=K0 1 1 

Uh - —Ta \th + 

- 

f(^) I 
1 : 

Kq ' 




-KO 


(23) 


Fig. 1^ shows how the cancellation between \th and \th breaks down for mo 7 ^ 0. We 
see in the left panel that the triangle and box fnnctions do not cancel when the scalar has a 
soft mass term, and the total amplitnde at threshold tends to a non-zero valne in the limit of 
inhnitely heavy scalar mass. The right panel shows how sensitive the cancellation between 
Fa'^ \th and F^'^ \th is to the presence of soft scalar mass terms, for a selection of fixed 
physical scalar masses. In addition to the cancellation ai k = ko, the amplitnde obviously 
vanishes when the scalar does not couple to the Higgs, k = 0. Discounting this trivial case, 
the amplitude quickly grows as we move away from the scenario where the scalar gets all of 
its mass from the Higgs. 
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FIG. 5; 2h partonic cross section in the SM, as a function of the Higgs self-coupling. Only top 
quarks are ineluded, with nit =173 Ge V. 

Finally, we examine the cancellation away from threshold. In addition to the spin- 0 
amplitudes realizing their full functional dependence on the partonic CM energy beyond 
s = 4m^, there are spin-2 contributions Gn to 2h production. The full amplitudes are 
known in terms of loop integrals, and reveal the strong cancellation near threshold when 
evaluated numerically. In the fermion case. Fig. shows how the partonic 2h cross section 
changes above threshold for the SM top. Near threshold, there is a pronounced cancellation 
between the triangle and box diagrams for the value of the Higgs self-coupling predicted by 


Eq. 22 This dip shifts and becomes much weaker as we move above threshold. 

It is interesting to consider how well our closed form expressions for the threshold 2h 
amplitudes approximate the full amplitudes. In the threshold amplitudes above, 

the Higgs mass may be considered as a proxy for the CM energy, s = 4m^, except in the 
triangle diagram where it appears in an s-channel propagator. This motivates the closed 
form approximation. 


F^\s,i,ml,m}) ^ \th (f 16mj/s) 
^ F^'^ \th (f ^ 16m,Vs), 


(24) 


for the 2h amplitude beyond threshold. Fig. shows the result of using this approximation 


in Eq. 15 for loop fermions and scalars of various masses. While the total cross section 


predicted by the expression in Eq. ‘M is not close to the true cross section, the normalized 
invariant mass distribution is fairly well reproduced. In particular, for rrihh > 2mj, two 
different loop particle propagators may go on shell, causing a nonzero imaginary piece of the 
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FIG. 6: Normalized invariant mass distributions for 2h production with a single loop particle. The 
upper panels show fermions of masses 173 GeV (left) and fOO GeV (right), while the lower panels 
show scalars of masses 173 GeV (left) and fOO GeV (right). The solid and dashed lines show the 


exact distributions and the approximation of Eq. 24 , respectively. 


exact amplitude that is visible as a feature iu the iuvariaut mass distributiou. Similarly, the 
threshold approximatiou iu Eq. 


24 


iucludes a term proportioual to arcsiu^ note 

that there is iu priuciple also a discoutiuuity iu the iuvariaut mass distributiou at rrihh > 4mj, 
captured by terms proportioual to arcsiu^ threshold approximatiou. However, 

there is uo visible correspoudiug feature iu the iuvariaut mass distributious of Fig. The 
discoutiuuities at mhh > 2 mi aud rrihh > dm* correspoud to the discoutiuuities iu the 
threshold amplitudes at r* = 4 aud r* = 1, respectively. Fig. shows the imagiuary parts 
of the threshold amplitudes, allowiug us to compare the discoutiuuities. The uuphysical 
discoutiuuity at = — 1 , which may be traced back to the arcsiuh terms iu the threshold 
amplitudes of Eqs. 17 aud 20, is uot showu. For both fermious [TU] aud scalars, we observe 
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FIG. 7: The imaginary parts of the 2h threshold amplitudes with loop fermions (solid) and scalars 
(dashed). Discontinuities occur when loop propagators go on shell. Appendix^ derives the imagi¬ 
nary part of the threshold amplitude with a loop scalar directly using cut techniques. 

a much larger discontinuity at r* = 4 than at r* = 1, due to the lack of any imaginary piece 
of the amplitude for r* > 4. The much smaller discontinuities in the imaginary part of the 
threshold amplitudes explain the lack of any visible features at rrihh = 4mj. 



IV. EXAMPLES 

The previous section demonstrated how the threshold cancellations between triangle and 
box diagrams render 2h production extremely sensitive to non-SM couplings. In this sec¬ 
tion we consider modihcations of the 2h distributions from anomalous fermionic Yukawa 
couplings, from colored scalar loops, and from fermionic top partners. Effects of anomalous 
fihh couplings on the kinematic distributions have been examined in [2H]. In some cases the 
allowed new interactions are severely restricted by the requirement that Ih production occur 
at the observed rate. In addition, we are interested in whether 2h production distributions 
can distinguish between fermion and scalar loop contributions. 

In all of our numerical results we use CT12 NLO PDFs [SOI IRT] with the associated NLO 
values for Og, and take = 173 GeV, rrih = 4.3 GeV, and rrih = 125 GeV. For Ih production 
we take n = hr = hf = rrih, while for 2h production we set n = hr = hf = rrihh- We use 
the LO 1-loop predictions for both Ih and 2h production. In addition, the 1-loop functions 
are evaluated using the software LoopTools |82], as well as independent in-house routines. 
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FIG. 8: Total cross sections for Ih and 2h production with anomalous top-Higgs couplings, nor¬ 
malized to the 1-loop SM prediction. The SM rate eorresponds to St = 0. All other couplings are 
assumed to be SM-like. 

A. Anomalous Yukawa Couplings 

We begin by considering the effects of anomalous top Yukawa couplings in Eq. [T} assuming 
all other couplings are SM-like. In Fig. we show the both the Ih and 2h rates, normalized 
to the one-loop SM rate as a function of the top quark Yukawa. For positive St, the require¬ 
ment that \Rh — 1| < .20 only allows a ~ 40% deviation in a{gg —)■ hh)/a{gg —)■ hh)sM-^ 
As measurements of the Ih rate become more precise, the allowed deviations for the 2 h 
rate due to an anomalous top Yukawa coupling will also become smaller. An important 
assumption throughout this work is that there are no light particles which could allow for 
resonant production of two Higgs bosons, in which case the rrihh spectrum would exhibit a 
clear peak at the mass of the resonance. The effect of a non-SM top Yukawa coupling on 
the invariant mass distribution is shown in Fig. i for both \SS = 13 and 100 TeV.^ When 
St % 0, the cancellation between box and triangle diagrams described in the previous section 
is spoiled and the resulting cross sections vary by up to a factor of two [S3l El]. This same 
variation is seen at both y/S = 13 TeV and 100 TeV. The effect of changing the top Yukawa 
coupling and the tri-linear Higgs coupling in a correlated manner can be quite dramatic, as 

^ We note that negative 6t is now excluded by global fits to Higgs couplings [HI |H] . 

There is not much difference in the invariant mass spectra between yfS = 13 and 100 TeV. This feature 
has also been observed in Ref. [25] . 
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FIG. 9: Invariant mass distributions for 2h production with anomalous top Yukawa coupling. The 
SM rate corresponds to 5t = 0. All other couplings are assumed to be SM-like. 


pp->hh, VS= 13 Te V pp^hh, Vs=100 Te V 
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FIG. 10: Invariant mass distributions for 2h production with anomalous top Yukawa coupling, 5t, 
and anomalous tri-linear Higgs couplings, 63 . The SM corresponds to St = 63 = 0. All other 
couplings are assumed to be SM-like. 


shown in Fig. Note the interesting cancellation for large and positive 63. 

In the SM, the contribntion of the b qnark is small for both the Ih and 2h prodnction [85] . 
although for large enough 5^ the production from b quark initial states becomes important 
[Sni EZ] • As the b quark Yukawa is increased, the rate for gg ^ h is substantially altered, 
while the 2h rate is rather insensitive to the b Yukawa as seen in Fig. Note that 

r{h —>■ bb) ~ (1 + 6byr{h —)■ bb)sM and the LHC experiments limit the total Higgs width 
l89] . Th^tot < so this implies a rough limit 6 b < 0(1). 
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FIG. 11: Total cross sections for Ih and 2h production with anomalous bottom Yukawa couplings, 
normalized to the LO SM prediction. The SM rate corresponds to Sb = 0. All other couplings are 
assumed to be SM-like. 

B. Fermionic Top Partners 


In this subsection and the following, we consider fermion and scalar contributions to 
gg —hh and pose the question: 

• Can we determine the nature of the loop particle, by examining the properties of the 
scattering amplitude? 


We use the analytic properties of the amplitude discussed in Sec. |^to draw some conclu¬ 
sions. 

We begin by considering the effects of a heavy color triplet fermion, with mass Mt, with 
a SM-like Yukawa coupling to the Higgs boson. The invariant mass distribution for this 
heavy fermion is compared with that from m* = 173 GeV in Fig. The distribution has 
an interesting dip near rrihh ~ 2Mt due to the presence of a cut in the amplitude. This dip 
persists even when the Higgs tri-linear coupling is allowed to have a non-SM value, ^3 7 ^ 0. 
At y/S = 100 TeV, the invariant mass spectrum has a more significant support at large rrihh, 
compared to that at y/S = 13 TeV. 

A 4th generation of chiral fermions would increase the rate for Ih production by roughly a 
factor of 9 above the SM prediction, far above the allowed region from current data [USl ES] • 
So we consider the addition of a vector-like fermonic top partner. In the simplest example, 
a top partner singlet model, there exists a charge-| SU{2)l singlet particle, which 
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FIG. 12: Invariant mass distribution for a heavy fermion with mass, Mt, and SM-like Yukawa 
couplings. The tri-linear Higgs coupling is allowed to vary from the SM value of 83 = 0. 

mixes with the SM-like top quark, T^. The Yukawa couplings in the top partner sector are 


— Ly 


rsj 


h/mn+hi-lun ++ i‘.c. , 


(25) 


where the Standard Model-like particles are denoted as 


fjL = 




n, hR 


(26) 


The addition of the A 5 Dirac fermion mass term in Eq. 25 means that the fermion masses 
are not completely determined by electroweak symmetry breaking. We can always rotate 
such that A 4 = 0 and so there are 3 independent parameters in the top sector, which we 
take to be the physical charge-| quark masses, rrit and Mt, along with the mixing angle, 9l. 
In the following, we will abbreviate sl = siiiOl, cl = cos9l. The couplings of the physical 
heavy charge-1 quarks to the Higgs boson are. 


mt 2 - , Mt 2 7Ti m , , Mt- mi. '^tm 


— Ch — —+ 

V V 


sfTiTiih + slCl - tiTjih + slCl—+ h.c.. (27) 

V V 


The parameters of the fermonic top partner model are limited by electroweak precision 
measurements to sin^t: < -12 and by direct search experiments to Mt > 880 GeV 

[931 EU. 
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FIG. 13: Ih and 2h production, normalized to the SM rate, in the top partner singlet model with a 
heavy top partner mass of Mt = 800 GeV. 

In the mt, Mt oo limit, single Higgs production in the top partner model is virtually 
identical to that in the SM ^ 


( 28 ) 


a 


/ 2 2 

top partner ^SM , 

= (^lh + t - 


Ih 


2’M2 


\mt .y.T. 

In Fig. we compare the rate for Ih and 2h production in the singlet top partner model 
with Mt = 800 GeV, as a function of the mixing angle, c^. For values of cl allowed by 
precision EW measurements, the Ih rate can be seen to be indistinguishable from that of the 
SM. On the other hand, 2h production receives contributions from the mixed tTh couplings 
of Eq. 1^ and in general can be quite different from the SM prediction, as shown in Fig. 
13 Even imposing the restrictions from precision EW data, 2h production can be reduced 
by up to 20% from the SM prediction |93l[9l], although the rate cannot be increased in this 
class of model. The relative reduction of the 2h rate is roughly the same at -sfS = 13 TeV 


and 100 TeV. The invariant mass spectrum are shown in Fig. 14 for the top partner model 


Because the EW precision constraints require that the mixing angle be very small, the rrihh 
distribution in indistinguishable from that of the SM. It would be interesting to investigate 
slightly less simple models by including tthh and TThh couplings, which arise in models 
where the Higgs arise as a pseudo-Nambu-Goldstone boson 


This is simply a statement of the decoupling limit. 
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pp^hh, VS=13 TeV pp^hh, Vs=100 TeV 

CT12PDFS 1^=1"^’ CT12PDFs 



FIG. 14: Invariant mass distribution in the SM (solid) and in the top partner model with = .9858 
and Mt = 800 GeV. 

C. Scalar Top Partners 


In this section, we compare the results of the previous section where the loop particles 
are fermions with rrihh distributions where instead of a fermion, there is a colored scalar in 
the loops. We begin by replacing the SM top quark in the gg ^ h triangle diagram and 
in the gg —)■ hh triangle and box diagrams with a color triplet scalar of the same mass. 


rris = mt = 173 GeV. Fig. shows the ratio of the total cross sections for both Ih and 2h 
production, normalized to the lowest order SM predictions, in this scenario. In the case of 
a color triplet scalar of mass = 173 GeV in the loops, we see that, in order to reproduce 
the SM rate for Ih production (the black dashed line), n need to be quite large, n < 2. If k 
is tuned to obtain a/asu = 1 for gg — h, then a color octet intermediate particle replacing 
the top quark with positive k (the solid black line) would predict a highly suppressed rate for 
2 h production (the red dashed line). Alternatively, we can tune both k and the scalar mass 
such that both Ih and 2h production have the SM rates, as shown in Fig. Although the 
total rates are identical to the SM predictions, the kinematic distributions from color octet 
and triplet intermediate states are quite different than those from the SM top, as plotted in 


Fig. [T^ The scalar needs to be quite light to reproduce the SM rates, and the distribution 
is peaked at much lower mhh than the SM prediction. 

We also consider Higgs production in the presence of the SM top quark and a colored 
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,^^S =13 TeV 


fns=mtop 



K 


FIG. 15: Comparison of Ih (dashed black) and 2h production (blue dot-dash) to the SM rate, when 
the SM top quark is replaced by a color triplet scalar with mass, mg = 173 GeV. The solid black 
(red dashed) curves correspond to the ratios to the SM predictions for Ih and 2h with a color octet 
scalar replacing the top quark. 


scalar®. Assuming the top Yukawa is SM-like, an additional scalar receiving all of its mass 
from electroweak symmetry breaking would give an unacceptably large contribution to the Ih 
production cross section, regardless of its mass and SU{3) representation. This immediately 


follows from Eq. 10 A heavy color triplet scalar with k, = changes the Ih production 

rate by 54%. Lighter scalars and scalars in other color representations result in even larger 


deviations. Fig. shows the effects of color triplet and octet scalars on Ih production 
in the large mass limit, as functions of the proportion of the scalar mass coming from the 
Higgs held. Heavy scalars receiving all their masses from the Higgs have mo = 0, and are 
not compatible with a simple average of current ATLAS [65] and CMS [66] limits on the Ih 
production rate from gluon fusion, which is drawn as a shaded band in Fig. [T^ 

The limits imposed from Ih production constrains the sensitivity of 2h measurements to 
reveal new scalars, especially those with masses close to the weak scale. This is because a 
light scalar, in addition to the SM top, modihes the Ih rate signihcantly, unless its coupling 
to the Higgs is small, which at the same time diminishes its impact in 2h kinematic distri¬ 
butions. However, heavy scalars will decouple quickly in the Ih rate and may show up in 


® In supersymmetry there are two colored scalars, the top squarks, mediating in the loop. Such a possibility 
is beyond the scope of current work and will be pursued elsewhere [98]. 
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8:m5=101.3 GeV; 3:m5=90.2 GeV; 'fs =13 TeV 



FIG. 16: Ratio of Ih production to the SM prediction when the top quark is replaced by a color 
triplet scalar of mass mg = 90.2 GeV (black dash) and by a color octet scalar of mass, = 101.3 
GeV (black solid), compared with the ratio of 2h production to the SM prediction when the top 
quark is replaced by a color triplet scalar of mass mg = 90.2 GeV (blue dot- dash) and by a color 
octet scalar of mass, mg = 101.3 GeV (red dash). 


iT{pp^hh)=crsjyi; Vs =13 TeV 



FIG. 17: Distributions for 2h production when the parameters are tuned to give the SM total cross 
sections for Ih and 2 h production. 


the high rrihh tail of the 2h distribution. The invariant mass distributions for 2h production 


are shown in Fig. 19 assuming a SM-like top quark and an additional 800 GeV color triplet 
scalar. If the scalar receives half of its mass squared from electroweak symmetry breaking, 
m-o = the Ih rate is in roughly 2 a tension with the current measurement, and the 

2h distribution deviates from the SM expectation starting at 2ms, roughly speaking. For 
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PP^h 



FIG. 18: Ratio of Ih production to the SM prediction with a colored triplet (blue solid) or octet 
(red dashed) scalar in addition to the SM top quark. The shaded band is the 2a limit on the rate 
obtained from an average of current LHC results. 


pp^hh.-vG" = 13TeV 
fi = mhh, CT12 PDFs 


pp^hh, VT = 100 TeV 
= ffihh- CT12 PDFs 



FIG. 19: Invariant mass distribution in 2h production with the SM top quark in addition to an 800 
GeV color triplet scalar that gets all (red dashed) or half (blue dot-dashed) of its mass from the 
Higgs. The SM (black solid) is shown for comparison. 

comparison, if the entire mass of the scalar was due to the Higgs, the feature at nihh = 2ms 
would be quite signihcant. 


V. CONCLUSION 

The observation of double Higgs production will be an important milestone in under¬ 
standing the nature of electroweak symmetry breaking. In the past the focus of this channel 
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has been on extracting the Higgs tri-linear self-coupling. In this work we showed that the 
goal can be much broader and encompass understanding the nature of the UV physics giving 
rise to Higgs coupling to two gluons, which is otherwise difficult to probe in single Higgs 
production. 


In Section |IH[ we examined the differences in the threshold behavior of double Higgs 
production resulting from intermediate scalar and fermion loops and in Section |IV[ we 
demonstrated that even if the parameters in a model with colored scalars are tuned to 
reproduce the SM rates for single and double Higgs production, the resulting invariant mass 
distributions can be signihcantly different from the SM. These distributions are also very 
sensitive to whether an additional scalar gets all of its mass from electroweak symmetry 
breaking. While Higgs plus jet production is also sensitive to the spin of loop particles and 
has a greater cross section, it does not enjoy the same large amplitude cancellation present 
in double Higgs production [72]. We also investigated the effects of anomalous top and 
bottom Yukawa couplings and showed that the resulting changes in single and double Higgs 
production relative to the SM rates are roughly the same at y/S = 13 and 100 TeV. 

Clearly, it will be an important experimental question on how to extract the wealth of 
information contained in the double Higgs production. Our work provides strong motivation 
to pursue this issue experimentally. 
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Appendix A: 2h Amplitudes from scalars 


Here we collect the contributions from virtual scalars computed in Ref. 
2 r 


jpis) _ _2^ 


— —A2Z1. 
^ ~ S 


.ml 


1 + 2m^C'o(s,m^) 


i, u, ml) + Do{s, it, i, ml) + Do{i, s, it, ml) 


(A.l) 


mi-u „ 2. ml-i . 2n , Pt 


= -4^<j-Go(s,m^,m2) 


-Go(M,m;) + 


^ Go(t, mj + s, M, mj 


(A.2) 


+ml ( Do{s, i, u, ml) + Do{s, u, i, ml) + Do{i, s, u, ml) 

' —2u{u — ml)Co{u, ml) — 2i{i — ml)Co{i, ml) + s(s — 2ml)Co{s, ml) 

(A.3) 


2{ut - ml) _ 

+su^Do{s, u, i, ml) + sPDq{s, i, u, ml) + s{s — 4m^)Go(s, ml, ml) 


where = {itt — ml)/s. In the above we have 


Cois,ml) = 


C{0,0,s,ml) 
f dPk 1 


J [A;2 _ ^2) + p^)^ — ml] 

Co{s,ml,ml) = C{ml,ml,s,ml) 

_ r d^k 1 1 1 

J ZTT^ [/c^ — ml] [{k + — ml] [{k + ki + ^ 2 )^ — ml] 

Do{s,i,u,ml) = D{0,0,ml,ml,s,u,ml) 

r d^k 1 1 1 


J ZTT^ [/c 2 _ ml] [{k + Pi)2 - ml] [{k + pi + ^ 2 )^ - ml] [{k + Pi + P 2 - - ml] 

Do{s,u,i,ml) = D{0,0,ml,ml,s,i,ml) 

_ r d^k 1 1 1 1 

J iTi^ [k‘^ - ml] [{k + pi)2 - ml] [{k + pi + ^ 2 )^ - ml] [{k + pi + P 2 - ^ 2 )^ - ml] 

Do{t,s,u,ml) = D{0,ml,0,ml,i,u,ml) 

_ r d^k 1 1 1 1 

J ZTT^ [/c^ — ml] [{k + pi)2 — ml] [{k + pi — /ci)^ — ml] [{k + pi + P 2 — ki)^ — ml] 

(A.4) 


and (pi + p 2 )^ = (^1 + ^ 2 )^ = s, {Pi - ki)‘^ = i, (Pi - ^ 2 )^ = u, and kf = kj = m 
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FIG. 20: Diagrams contributing to Fr- Bubble diagrams with quartic scalar-gluon vertices are 
included. The cuts shown are used to calculate 


Appendix B: Closed Form Amplitudes for gg 


hh 


Here, we calculate the imaginary part of the amplitude for 2h production from scalar loops 
at threshold, using cut techniques. This is a new result analogous to the recent computation 
for fermions |73|. By using the dispersion relation, we may recover the full 2h amplitude in 


closed form, which is then analyzed in Section III 


We start with the general amplitude in Eq. For simplicity, we assume a single scalar 
that gets all of its mass from the Higgs, so that mo = 0 and k = kq = 2m^/u^ in Eq. 
Our results can easily be generalized to scalars with arbitrary couplings and masses, and 
we emphasize that they do not assume a heavy loop particle. At threshold, s = 4m^, only 
the spin 0 piece contributes [61]. The imaginary parts of the corresponding form factors 
can be obtained from cutting all possible gg —)■ hh diagrams, and sending all cut 

(s) ( 5 ) 

propagators on shell. We will compute ImF^'^ and ImT^ ^ separately at threshold. 

I I (s) 

Fig. 20 shows the diagrams that are responsible for the form factor. In addition to 


the triangle diagram which may be obtained by replacing the top quark in the SM double 
Higgs triangle diagram with a scalar, we have the additional s*sh‘^ coupling. We also include 
bubble diagrams with quartic scalar-gluon couplings with the above diagrams in the triangle 
form factors, as they are related through gauge invariance. Now, the imaginary part of the 
double Higgs amplitude receives contributions from the cuts shown in the diagrams of Fig. 
20 through 

\mM D j dU2MLMR (B.l) 

where A4, A4l and A4r refer to the full double Higgs amplitude and the left/right halves 
of a cut diagram. The integral f dl \.2 is over the phase space of the cut propagators. Each 
cut diagram in Fig. 20 contributes separately to The halves of the cut diagrams 


are simply tree-level amplitudes for gg —>■ s*s and s*s —?• hh. Furthermore, since we are 
interested in the amplitude cancellation at threshold, we may project out the spin 0 piece 
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FIG. 21: Diagrams contributing to and Gg . Triangle diagrams with quartic scalar-gluon 
vertices are included. The cuts shown are used to calculate ■ 

of the amplitude to get \th- The kinematics of Appendix simplify considerably for 

s = 4m^, and we are left with 

1 + \/l ~ g/4 


Im F. 


W I _ tt / 4 


^ U - d ‘"“A 


v'l -tJA 


Cut diagrams that contribute to and are shown in Fig. 


21 


(B.2) 


The first two 


diagrams of Fig. 21 have identical cuts to the diagrams of Fig. 20, and the left sides are the 
same as in the earlier diagrams. At threshold with the cut propagators on shell, comparison 
of the right sides of these diagrams with those of Fig. immediately gives 


ImF^'^ \thD -^ImF^^ \th 


This is the contribution of the top row of Fig. 


21 


to IraF^^ \th. 


The contributions of the cuts in the third and fourth diagrams of Fig. 


21 


(B.3) 


to ImFn'^^ at 


threshold may be computed from the tree-level amplitudes for gg —)■ s*sh and s*s —?• h. 
Note that the two adjacent propagators attaching to either external Higgs may be cut, each 
choice leading to an identical set of contributions to the imaginary amplitude. Only one 
such set of cuts is shown in these diagrams. Both sets of cuts together yield 




a/T 


i + OiogiiHZl 
2I 


(B.4) 


Also, there is no contribution to the imaginary part of the gg —)■ hh amplitude from cutting 
two adjacent propagators attaching to an external gluon, because the amplitude for g ^ s*s 
is zero when the scalars are put on shell. 

Finally, the last cut diagram of Fig. 


21 


gives a contribution to Im F^ ^ that may be 


calculated at threshold from the gh —)■ s*s amplitude. We proceed as before, and find the 
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final contribution 




□ ' \th^ ^T{s)e - l] t‘ 1 {\/l + Ts- (l + log ^ 


16 


To 


1 — a/I + Ts 


(B.5) 


The sum of the right-hand sides of Eqs. B.2 B.3 , B.4 and B.5 gives the full imaginary 
2h amplitude at threshold. Now, we turn to the limits of the full amplitude as —)■ 0, oo. 
In the limit Tg —)■ 0, the amplitude vanishes since for a scalar that gets all its mass from the 
Higgs, K is proportional to through 


K = 


rr^ 

2v^ 


Tg 


(B.6) 


On the other hand, in the inhnite scalar mass limit —)■ cxd we may apply the low-energy 
theorem. From the effective Lagrangian for the interaction between scalars and gluons m, 
we know that the 2h amplitude goes as 


Ahh oc ( hh 


= ( hh 


2h h^ 
log 1 + — + ^ 

/J1 /JlZr 


0 


2h 

V 


/ 1.2 


0 


(B.7) 


which vanishes for the SM Higgs self-coupling | 721 . Given the limiting behavior of the 
amplitude combined with full knowledge of its imaginary part, then, the dispersion relation 


gives the full amplitude in Eq. 20 
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